
setwd("/Users/italolopez/Dropbox/NEP Evaluation - The role of beliefs/JPE Final Version/JPE Submission/Replication package/Outputs/Tables")
############################################################################
############################################################################
### DV: zptevir_irtscore2 ###
############################################################################
############################################################################

################# MV: Behavioral Index ###################
set.seed(7688);
require(MASS)
a=0.108
b=0.048
acov <- matrix(c(
  0.002116, 0,
  0,0.0004
),2,2)
rep=20000
conf=95
pest=c(a,b)
mcmc <- mvrnorm(rep,pest,acov,empirical=FALSE)
ab <- mcmc[,1]*mcmc[,2]
low=(1-conf/100)/2
upp=((1-conf/100)/2)+(conf/100)
LL=quantile(ab,low)
UL=quantile(ab,upp)
LL_perc=(LL/0.099)*100
UL_perc=(UL/0.099)*100
LL4=format(LL_perc,digits=4)
UL4=format(UL_perc,digits=4)
print(paste("T9-CIs","LB%",LL4,"UB%",UL4 ,sep=";" ))
out <-c("T9-CIs",LL4,UL4)
capture.output(out, file = "T9-CIs.txt", append = FALSE) 
M=mean(ab)
med=median(ab)
